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Statistical mechanics of confined binary system: 
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ABSTRACT 

We have used a toy model to study the behaviour of confined binary systems in 2D and 
O^ ■ compare it with previously known result in 3D. In the case of 2D, in which canonical 

0^ I distribution exists only above a critical temperature, we evaluate the exact form of 

partition function for this system and compare the exact partition function with the 
CJ I mean field partition function for the case of two-particle system. In contrast of its 

^ . 3D counterpart, there is no phase transition here. If this system, however, studied in 

^H I microcanonical ensemble, it shows two different phases of kinetic energy dominated 

f— ^ I with positive specific heat and potential energy dominated with negative specific heat 

Cr^ I in presence of short distance cutoff. In absence of short distance cutoff, surprisingly, 

the negative specific heat region will be replaced by region of large specific heat. This 

feature is completely new and there is no such a case in 3D. 
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1 INTRODUCTION 



O ■ The statistical behaviour of N particles interacting through Newtonian gravitational forces is very different from the statistical 
' [ behaviour of other many body systems such as neutral gases and plasmas. The central feature of gravitating system, in contrast 
» I • to normal many body systems, is the unshielded, long range nature of gravitational force. Because of this feature one of the 
Ti fundamental concepts of statistical mechanics, the extensive nature of energy, breaks down. This, in turn, leads to different 
C^ physical descriptions for the gravitating systems in the microcanonical and canonical distribution. 

^ _ The statistical behaviour also strongly depends on the spatial dimension. For instance, in 3D, the available phase volume 

for the system diverges and one is forced to use short distance cutoff. However the situation in 2D is different. In this case 
there is a microcanonical description for all values of energies, through the canonical approach exists only above some critical 
fS ■ temperature. 

We shall study some properties of these gravitating systems by introducing a toy-model (originally used in Padmanabhan, 
1990), based on a simple Hamiltonian, describing two particles of finite size, confined inside a box. This system shows several 



important properties of more complicated systems studied for example in D. Lynden-Bell and R. M. Lyndel-Bell (1977). We 
will study this toy model in both the microcanonical and canonical approach. 

In section 2 we will overview the properties of the toy model in 3D, [which was earlier done in Padmanabhan, (1990)] 
for providing the background needed for comparison with the 2D case, which will be studied in section 3. In section 4, we 
compare the nature of this system in 3D and 2D. 

An interesting feature of this 2D system is that the thermodynamic functions are all calculable analytically. This contrasts 
with the thermodynamical behaviour of 3D confined binary system. By studying the "isothermal cylinders" we find that these 
systems are remarkably similar to the simple toy model and that the system cannot exist at T < Tc where Tc is given by 
Tc = (l/2)Gm} . We also find that by introducing a short distance cutoff, in contrast to 3D case, the specific heat of the 
system would become negative in some intermediate temperatures. 
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2 OVERVIEW OF 3D CASE 



In this section we introduce and review a toy model which was first given by Padmanabhan (1990). In this model, by 
constructing the "statistical mechanics" of two particles of finite size confined inside a box and interacting via 3D gravity 
are described. This system exhibits several important properties of more complicated gravitating systems in spite of the fact 
that it has only two particles. In particular, this system exhibits the following two features, which seem to be generic to all 
gravitating systems: (i) When studied using the microcanonical ensemble, the system shows evidence for two different phases: 
a high temperature phase, dominated by kinetic energy, and a low temperature phase dominated by the potential energy and 
stabilised by some short distance cutoff which is of non-gravitational origin. Both these phases have positive specific heat. 
These two phases are connected at intermediate temperatures by a region of negative specific heat; this is precisely the range 
in which the kinetic and potential energies of the system are comparable and the system is in virial equilibrium, (ii) If the 
same system is studied using the canonical ensemble, the intermediate region of negative specific heat is replaced by a sharp 
phase transition. 

Then, in order to show that this toy model exhibits several important properties of more complicated gravitating systems, 
we shall compare the equilibrium of a gravitating system in the mean field limit, by evaluating the partition function for this 
system. In the absence of any short distance cutoff to gravitational interaction, the mean field solution is given by an isothermal 
sphere. It turns out that this system is remarkably similar to our simple toy model studied earlier. In fact, this mean field 
analysis confirms all the conjectures based on the toy model. 

2.1 Microcanonical description 

We will begin by studying the statistical mechanics of our toy model described by the Hamiltonian 

where (Q, P) are coordinates and momenta of the centre of mass, (r, p) are the relative coordinates and momenta, M = 2m is 
the total mass, /j, — (ra/2) is the reduced mass and m is the mass of individual particles. The range of r varies in the interval 
(a,R). This is equivalent to assuming that the particles are hard spheres of radius (a/2) and that the system is confined to a 
spherical box of radius R. In microcanonical distribution corresponding to this toy model (which is the relevant ensemble to 
the gravitating systems ) we define the entropy S[E) and temperature T{E) through the relations 

S{E) = \ng{E), T{E) = {dS{E)/dEr^ = {d\ng{E)/dE)-^ . (2) 

For our system the phase volume g{E) of the constant energy surface H — E is 

^R^i-Ey^il + aE/Gn?f, -Grn^/a < E < -Gm'^/R, 

9iE) _} (3) 



' iR^i-Ey^Kl + RE/Gni^f -{1 + aE/Gm'^f], ~Gm^/R<E<oo 

This function g{E) is continuous and smooth sA, E — {~Gm^ /R). 

Now the thermodynamic properties of the system can be analyzed from the T{E) curve. 

In the range when —Gm^ /a < E < —Gm^ /R we can write T{E) in the dimensionless form 

where t = {aT/Gn?) and e = [aE/Grn^). 

At e = — 1 which corresponds to the lowest energy admissible for the system, the dimensionless temperature t vanishes. 
It is obvious that in this case of e ~ — 1, i(e) dominated by the first term of (^j. As we increase the energy of the system, the 
temperature increases, which is the normal behaviour for the system. This trend continues up to 

£ = ei = -i(V3-l)~-0.36, (5) 

at which point the t{e) curve reaches a maximum and turns around. As we increase the energy further the temperature 
decreases. The system exhibits negative specific heat in this range. 

As one can see from (W), for realistic systems, i.e., R^ a only a small region in the range of —{Gm? /a) to — (0.36Gm^/a) 
we will have positive specific heat; for the rest of the region the specific heat is negative. In fact the existence of the positive 
specific heat region is due to nonzero short distance cutoff. In the absence of this nonzero short distance cutoff, the first term 
in (0) will vanish and we will get t ex — e"^ and negative specific heat in this entire region. 

For high energies limits, i.e., E > —Gm^/R, the second expression in (13) for g{E) will give 

'3{il + e)^^{R/a)[l + iR/a)e]^} l""' 



'^^^^l (l + .)3-[l + (i?/a)£]3 S) ■ (6) 
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E = (Ea/Gm^) 

Figure 1. The function T(E) for a confined binary system. 



This function will match with (Uh at £ = ~{a/R). It will decreases as we increase the energy, for a while, and then very soon 
it starts to increases as energy increases at some e = £2- Thus system will enter another phase with positive specific heat. The 
form of i(e) is shown in figjiJ. The specific heat is positive along the portions AB and CD and is negative along BC. 

For E ^ E2 = —{Gim?/R), gravity is not strong enough to keep r < R and the system behaves like a gas confined by 
the container; we have high temperature phase with positive specific heat. As the energy decreases to i? < E2, the effects of 
gravity begin to be felt. For E\ = —(Grn^/R) < E < E2, the system is unaffected by either the box or the short distance 
cutoff; this is the domain dominated entirely by gravity and we have negative specific heat. As system goes to E ~ Ei, the 
hard core nature of the particles begins to be felt and gravity is again resisted. This give rise to a low temperature phase with 
positive specific heat. 

It is also interesting to study (i) the effect of increasing R, keeping a and E fixed, and (ii) the effect of varying a, keeping 
R fixed. In former case, it is amusing to note that, if 2 < R/a < (v^ + 1), there is no region of negative specific heat. As we 
increase R, this negative specific heat region appears. In the latter case, at first one should rescale variables using (Gm^/R). 
This can be easily done and fig.0 shows the behaviour of the T{E) curve as the lower cutoff a is changed. As a is lowered, 
the negative specific heat region becomes more and more pronounced. If a is zero, we have negative specific heat for all 
E < -{Gm^/R) (see fig.|). 



2.2 Canonical description 

It is of interest to look at our system from the point of view of the canonical distribution. To do this we have to compute the 
partition function 



Z{f3) = / d^P d% d'^Q d'^r exp{-f3H), 



(7) 



Z{t) = t^{R/af I dx x^ exp{l/xt), 



which after integrating over P, p and Q and omitting an overall constant, which is unimportant, in dimensionless form becomes 

(8) 

1 

where t is the dimensionless temperature defined before. One can show that Z{t) can be well approximated by the expression 
© 0000 RAS, MNRAS 000, 000-000 
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LABEL = Log,Q(R/a) 




Figure 2. The T{E) curve on the short distance cutoflf a. For a > (3^'^ + 1) ^R there is no region of negative specific heat. As a is 
lowered, the negative specific heat region becomes more and more pronounced. 




E = (ER/Gm ) 

Figure 3. The T{E) curve for a = 0. The low temperature region with positive specific heat does not exist in the absence of the short 
distance cutoff. The system also exhibits a lower bound on the temperature. 
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Log,o(R/a) = 2.5 




Figure 4. Comparison of the T{E) relations for the canonical and microcanonical distributions. The negative temperature region of the 
microcanonical distribution is replaced by a phase transition in the canonical distribution. The microcanonical temperature is replaced 
by a factor 2/3 for convenience of comparison. 



{R/aft'^{l~2t)~^ exp{l/t), ior t < tc 



(9) 



U^{R/af{l + 3a/2Rt), 



for t > tc 



where tc = [3ln{R/a)]^ , is the critical temperature at which the transition occurs. Given Z{f3) one can compute E{l3) by 
the relation -E(/3) = — {din Z/df3). This relation can converted to give T{E), which can be compared with the T{E) obtained 
earlier using the microcanonical distribution. From M we get 



aE/Gm? = 4t - 1, for t < U, 



(10) 



3i - 3a/2i?, 



for t > tc 



Near t ~ tc is a rapid variation of the energy and we cannot use either asymptotic form. The system undergoes a phase 
transition at f = tc absorbing a large amount of energy, 
1 



The specific heat is, of course, positive throughout the range. This is to be expected because the canonical description cannot 
lead to negative specific heats. 

The exact T(E) curves obtained from the canonical and microcanonical distributions are shown in fig.0. The descriptions 
match very well in the regions of positive specific heat. The negative specific heat region of the microcanonical distribution is 
replaced by a phase transition in the canonical description. 

The physics of canonical distribution is best understood by studying _B as a function of T. As we increase the temperature 
from zero to, the energy increases from the ground state value —{Gm? /a) in accordance with first equation of (ILCJ). As the 
temperature approaches Tc and cross it, a phase transition occurs in the system and the energy increases rapidly. The latent 
heat in the system is large enough to push the system into the high temperature phase. At still higher temperatures, the 
energy increases steadily with the temperature in accordance with second equation of (|lO|) . 

We can now compare the canonical and microcanonical descriptions for our system. At both very low and very high 
temperatures, the descriptions match. The crucial difference occurs at the intermediate energies and temperatures. The 
microcanonical description predicts a negative specific heat and a reasonably slow variation of energy with temperature. The 



© 0000 RAS, MNRAS 000, 000-000 



6 Ali Nayeri 




Figure 5. The u — v curve fore the isothermal sphere without any short distance cutoff. 



canonical description , on the other hand, predicts a phase transition with a rapid variation of energy with temperature. Such 
phase transitions are accompanied by large fluctuations in the energy, which is the main reason for disagreement between the 
two descriptions. 

In order to confirm all the conjectures made earlier based on the toy model we shall briefly explain the equilibrium of a 
gravitating system in the mean field limit. In the absence of any short distance cutoff to gravitational interaction, the mean 
field solution is given by an isothermal sphere. In this case we may write, 



V = AnGpce 



-/3[0(a:)-0d 



where pc = p(0). 

One may introduce the following dimensionless variables in order to write the isothermal equation, 

x = r/Lo, n = p/pc, m = M{r)/Mo, y = f3[(t> - (l>c]. 

where Lq and Mo are some characteristic length and mass, respectively. 
Then in terms oi y{x) the isothermal equation becomes 

x'^ dx \ dx J 

with the boundary condition y(0) = y'{0) = 0. By defining the new variables as 



m/x, 



3 / 2 / 

/m — nx fv, 



one can express the isothermal equation as the following: 

u dv M — 1 

V du u + V — 3 

with the boundary conditions t; = at u = 3 and dv/du — —(5/3) at (u, v) — (3, 0). 

The solution curve starts at (3,0) and spirals indefinitely around the point (1,2) as x tends to infinity (see fig.t 
isothermal spheres must necessarily lie on the this curve . 



(12) 



(13) 



(14) 



(15) 



(16) 



All 



The u — V curve also implies a bound on the temperature of the system [Lynden-Bell and Wood (1968)]. It is clear that 
for any isothermal sphere v is bounded from above , i.e, v < Vmax, where Vmax ~ 2.5. Since 



V = m/x = {M/Mo)iR/Lo) = {GM/R)f3, 



(17) 
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0.2 

(RE/GM^) 



Figure 6. The T{E) curve for the isothermal sphere without any short distance cutoff. The specific heat is positive along AB and 
negative along BC. 



we immediately get 
T '> T T 



qagm'^/r. 



(18) 



The T{E) curve for the isothermal sphere can be determined by combining the relations 



r = {GM''/R)v-\ 



E = {GM^/R)v-''{u- -). 



(19) 

The form of the curve is shown in figM The specific heat is positive along the AB (which represents a high temperature 
phase) and is negative along BC. The branch CD is unstable and is not physically realisable. Since we have not introduced 
any short distance cutoff in the system we do not have a low temperature phase with positive specific heat. If we do that, 
a short distance cutoff makes the curve turn to the left on the CD branch and produces a second region of positive specific 
heat. Figure ^ is analogous to fig.|3[ which describes the T{E) curve for a binary model in the absence of a short distance 
cutoff. The modification of fig.g in the presence of a short distance cutoff will be very similar to the modification of fig.^ into 
figJd (For example see fig. 4. 11 in Padmanabhan, 1990). 

3 THE CASE OF 2D BINARY 

We shall now derive corresponding results for the 2D confined binary system. We will see that the thermodynamic functions 
are all analytically calculable for this case. It is also interesting to contrast the thermodynamical behaviour of the 37? confined 
binary system with the different, but eventful, 2D one. 



3.1 Microcanonical approach 

Let us consider two particles interacting via 2D gravitational force and moving in a 2D region of radius R. The potential in 
2D gravity satisfies the Poisson equation: 

V'^cf) = 2TiGp. (20) 

For point particle <j) will be logarithmic; the potential energy of interaction between two such particles will be, 

(21) 



(/(xi,X2) = Grn^ In 



R 
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E = (ER/Gm^ Log(R/a)) 



Figure 7. The T{E) curve for the low energies region for 2D. In spite of existence of short distance cutoff, there is a region which 
specific heat is negative, although very small. 



We shall first study the microcanonical distribution corresponding to our system. The phase volume g{E), which is the 
volume of the constant energy surface, in the phase space: 



9{E)= lf[d\id'p^S{E^H), 



(22) 



where, H, is the Hamiltonian for a 2D gravitating system of 2 particles with logarithmic potential, defined as. 



where, as in the 3D case, (Q,P) are the coordinates and momenta of the centre of mass, (r, p) are the relative coordinates 
and momenta, M = 2m, is the total mass, fi — {m/2) is the reduced mass and m is the mass of the individual particles, we 
shall further restrict the range of the coordinate r to the interval {a,R). Where, a, is the short distance cutoff equivalent to 
assuming that particles are hard spheres of radius (a/2), while, R, is the large distance cutoff equal to the confining radius of 
the system. Therefore, the phase volume g{E) becomes. 



giE) = AR^ 



E - Gm^ In 



ii) 



dr. 



(24) 



With the range of integration in (p4|), limited to the region in which the expression in parentheses is positive, i.e., [E 
Grn^ In {r/R)\ > 0. Therefore, 



R exp [E/Gm^), -Gm^ \n{R/a) < E < 0, 



(25) 



R, 



< E <oo. 



A is some constant which is irrelevant to our discussion. 
The integration, then, will yield the following result: 



Gm? 



\AR^{R^e 



2 2E/Gm^ 



a^{2E/Gn? + 2 \n{R/a) + 1)}, 



-Grr? \n{R/a) < E < 0, 



(26) 



^AR^{R'^{2E/Gm'^ + 1) - a^{2E/Gm'^ + 21n(i?/a) + 1)}. Q<E <oo. 
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E = (ER/Gm^ Log(R/a)) 

Figure 8. The T(E) curve. Note that the negative specific heat region is not visible due to its sinallness compare to high energy region. 




c = (ER/Gm^ Eog(R/a)) 

Figure 9. The T(E) curve for difi'erent short distance cutoflf a. As one decreases R, the negative specific heat will become more and 
more pronounce. 
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E = (ER/Gm ) 

Figure 10. The T{E) curve for the special case, a = 0. The system does not exhibit any negative specific heat. Compare this figure 
with the fig. 3. 

It is obvious that g{E) is continuous and smooth at -E = 0. We shall study the thermodynamics of the system, which can be 
now studied using g{E), in the two regimes given above. 

The entropy of the system, S{E) — \ng{E) in the case of very low energies, i.e., Gm? \n{a/R) < E < will be. 



S{E) = In 



-AGm'^R^ 
4 



+ ln 



2 2E/Gm^ 



R e 



2E 



2 In 



+ 1 



and the temperature of the system, T{E)~ 



( 

dS{E)/dE with the help of (|27[) is. 



(!) 



T{E) = Gm^ 



[E/Gm^ + ln{R/a) 



or in dimensionless form as, 

1 (e+1) 



t(e) = 



+ 



c)^ 



(27) 



(28) 



(29) 



21n(i?/a) l-(i?/a)2(i- 

where we have defined f(e) — T{E)/Gm? \n{R/a), and e — EjGrr? ln{R/a). At the lowest energy admissible for our system, 
which corresponds to e = —1, the temperature is f(e) = [21n(_R/a)]~^. It is clear that in (|23) for e ~ —1 the first term 
dominates. So as we increase the energy of the system, the temperature decreases. This behaviour continues up to e = at 



which the point f(e) curve reaches to it's minimum, f(e) 
specific heat region for — 1 < e < (see figj?]). 

For _B > we should use the second expression in (k( 



l/{21n(_R/a)} — a /(R — a ). Therefore, we obtain a negative 



SiE) = In 



-AGm'^R^ 
4 



+ ln 



R- 



/ 2E 
[Gm? 



+ 1 



( 



2E 



\Gm' 



for g{E). In this case the proper expression for S{E) would be 
'R 



+ 2 In 



(!) 



and thus we get 

1 

t{e) = £ + 



2ln{R/a) 1- {R/a)^ 



(30) 



(31) 



This function, clearly, matches with (g9|) at e = 0. As we increase the energy , the temperature continues to increase. Thus 
in high temperature phase we enter to the positive specific heat region. The form of t{e) is shown in figjs] (Due to smallness 
of negative specific heat region compared to high energy limit, the negative specific heat region is not visible in the figure). 
However, as we decrease the long distance cutoff R, the negative specific heat region becomes more and more pronounce, see 
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figi 



If we set a — 0, we easily get, 



T{E) 



1/2, 



-cxj < £■ < 0, 



(32) 

1/2 + E/Gm'^, 0<E<oo. 

It is obvious that in the absence of the short distance cutoff, there will be no negative specific heat region ( fig .FloI). 
Comparing this result with it's counterpart in 3-D case for a = shows that negative specific heat replaced by constant 
temperatures. In fact, the effect of short distance cutoff in 2D case is just destablising the effect of gravitational potential 
energy, unlike its effect in 3-D potential. By putting a short distance cutoff, we actually distort the phase space. In other 
words, having a short distance cutoff in the system is equivalent to removing some part of phase space. This distortion in 
phase space may cause many unexpected consequences. In 2D, for instance, removing some part of phase space causes an 
unexpected region of negative specific heat, whereas for 3-D, we get region of positive specific heat. Comparison of 2D and 3-D 
with short distance cut off, though, suggests the E — T graph in 2D in many aspects is similar to its 3-D counterpart, except 
that in 2D there is no initial positive specific heat region. In other words, the E ~T graph in 2D is almost same as the E — T 
graph in 3D, but as if the graph has been shifted to the right: the starting point of the 2D graph almost coincides with the 
maximum of 3-D one. 



3.2 Canonical approach 

Let us now consider the partition function Z{/3) of the system, which is given by the integral, 

f +00 

-0E 



Z{I3) 



dE g{E)e~ 



/' + 00 



dE g{E)e'"'^ + 



dE g{E)e' 



Z1+Z2 



(33) 



The range of integration is from (—00) to (+00), since negative values of E are allowed. Thus in our case, that interval will 
break to two separate areas. We have for Zi: 



1 r° 

ZiiP) = -AR'^Gm'^ / dE exp ([2/Gm^ - f3]E) 

•/ — oc 

The above integrand will diverge in the lower limit, unless, 

Then we obtain, 

ZiW)--AR 2_^^2^. 

Since ^i(/3) diverges at Gw?Pc — 2, the system can exist only at /3 < f5c- As for Z2, 



(34) 



(35) 



(36) 



Z2W) 



\aR* I dE{Gm'' + 2E)e~^'^ 



- 2^U 



1 + 



/3Gm= 



(37) 



which is the Z{f3) as one can obtain from saddle-point approximation (Padmanabhan 1991). Thus, we see that the saddle-point 
approximation is accurate at high temperatures, i.e. for Gni^(3 ^ 1, and the mean field approximation must break down at 
low temperatures. If now we introduce the short distance cutoff, a, the corresponding Zi and Z2 will be modified as 

^ 1 . 

(38) 



Z,{P) = -AR 
and 



1 . „4 {Gm'f 



1 Ra 

2~Gm?(i 2 \ (5 



1 + -Gm^ [2 In (-R/a)-M]/? 



Z2(/3) = \aR* I dE{Gm'' + 2E)e-^'' 



- 2^U 



1 + 



f3Gm"^ 



I J Ra 
2^ - 



1 + -Gm'' [2 In (.R/a) + l]/? 



(39) 



Given Z{T), one can compute the mean energy of the system, which is given by 



E{T) =r {din Z/dT), 
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Figure 11. The T/Tc as a function of the total energy of two dimensional gravitating system. There is no region of negative specific 
heat but the system exhibits a lower bound on the temperature. 



or in dimensionless form as 



r(91nZ/ai), 



where for high energies, regardless of whether there is a short distance cutoff, the dimensionless energy scales as 



2r 



£2 



t. 



(41) 



(42) 



2t + l 

However at t ~ tc = l/[2\n{R/a)], energy diverges. Thus, again at Gm?Pc ~ 2, the partition function will blows up even in 
the presence of short distance cutoff a and in this case the system cannot exist at /3 > /3c, as well. In other words, the canonical 



description of the system exists only at sufficiently high temperatures ( Katz and Lynden-Bell 197q ). Note that there is no 
phase transition for this two dimensional system either with or without short distance cutoff. As the energy of the system is 
lowered, the temperature continuously decreases and asymptotically approaches Tc = (l//3c). 

Comparison of the canonical and microcanonical descriptions of our model shows that at very high temperatures, the 
descriptions match. The main difference occurs at low temperatures and energies. The microcanonical description predicts 
either infinite or negative specific heat, depend on whether there is a short distance cutoff, at low energies, whereas in canonical 
approach, there is no physical state below a critical temperature. Therefore, for low energies the two approaches disagree. 

As we compared the result of 3D binary with isothermal sphere, we shall outline the result of studying isothermal cylinders 
in order to confirm our earlier claim based on 2D binary. 

An isothermal self-gravitating cylinder is described by Poisson's equation, (EOl), in 2D, 



with p(r) given by 
p{r) = Aexp[~Pip{r)], 



(43) 



(44) 



where the constants A, and (3 are to be determined in terms of total mass M and energy E of the system. Regular solution 
of (Et3) have been given in Ostriker (1964), and Stodolkiewicz (1963): 



ip{r) = GM In ii+ 2/3" Mn 



1 



GM/3(1 - r^/R^) 



+ constant, 



(45) 



© 0000 RAS, MNRAS 000, 000-000 



confined binary system 13 

where R is the radius of the confining box. The potential on the axes, _R = 0, is real if (l/4)GM/3 is smaller than unity and 
becomes infinite in the limit (l/4)GM/3 = 1. Thus, there exist a lower bound, T^. Now all other physical variables like density 
p(r), pressure on the confining wall P{R) = I3~^p{r), and the energy 

E = M/3-1 + i / pipdx, (46) 



can be computed from (ktSh. We get 

E = jGM^ [2 In 7? + 2{T/Tc) + {T/T^f ln(l - Tc/T) + contstant] , (47) 

PV = N{T-Tc), (48) 

where Tc — {1/A)NGM'^ is the same critical temperature as that found earlier (for N = 2) in the canonical approach. It is 
clear the system can not exist for T < Tc, the pressure becomes negative, and the potential at the origin diverges. 

The T{E) curve for this system is shown in fig-O. It is clear that there is no region of negative specific heat for this 
two dimensional system. As E is lowered, the T{E) continuously decreases, however, T will never reach Tc since for that to 
happen, an infinite amount of energy has to be given away. As T tends to T^ the pressure on the wall tends to zero, see 
equation (|4^), the density in the centre grows continuously as 

p(0) = y(l-Tc/T)-\ (49) 

and the density contrast grows as, 

p{0)/p{R) = {l-TjTr\ (50) 

The system, after shrinking, collapses to thin dense string. Therefore, isothermal cylinders in contact with heat bath whose 



temperature T is slightly smaller than Tc are unstable and giving up an unlimited amount of energy ( Katz and Lynden-Bel] 
197q^ However, a comparison with the case of 3D binaries shows that, two-dimensional systems are more stable than three- 
dimensional ones. As we obtained in eq. (hsl) there also exists a Tc for isothermal spheres, in which there is no equilibrium. If 



an isothermal sphere has a density contrast less than 32 (Lynden-Bell and Wood 1968) and temperature slightly hotter than 



Tc, while the density contrast grows, it will lose energy to heat bath and cool down. As the density contrast keeps growing 



above 32, the specific heat becomes negative. When density contrast reaches the value 709 (Katz and Lynden-Bell 1978) and 
becomes very hot, the system is unstable and not physically realisable. Then after, isothermal starts collapsing in which the 
centre of the system becomes smaller and hotter whereas giving up energy to the outside parts of the isothermal sphere. 

Instabilities in isothermal spheres are mainly due to "wiggling" of potential around singular solution (see, e.g. fig. 6). There 
is no singular solution and thus wiggling of potential in 2D (see fig. 12). 



4 CONCLUSION 

It seems a binary system exhibits several important properties of more complicated gravitating systems in spite of the fact it 
has only two degrees of freedom. In particular, this system exhibits the following two features, which seems to be generic to all 
gravitating systems: (i) When studied using the microcanonical ensemble, the system shows evidence for two different phases: 
a high temperature phase, dominated by kinetic energy, and a low temperature phase dominated by the potential energy and 
(de) stabilised by some short distance cutoff in (2)3D which is of non-gravitational origin, mainly due to distortion of phase 
space. Both these phases have positive specific heat in 3-D whereas in 2D the latter phase has negligible negative specific heat 
for the case of {a/R) <C 1. In 3D case, these two phases are connected at intermediate temperatures by a region of negative 
specific heat; this is precisely the range in which the kinetic and potential energies are comparable and the system is in virial 
equilibrium, (ii) If the system is studied using canonical ensemble, the intermediate region of negative specific heat in 3D is 
replaced by a sharp phase transition releasing a large amount of latent heat. This suggests the following analogy: Gravitating 
systems in virial equilibrium are similar to normal systems (with short range forces) at the verge of phase transition. For the 
case of 2D, however, the canonical description leads to the completely different picture in low energies. The system does not 
exist below some critical picture Tc = (l/2)Gm^. 

On the other hand, the isothermal considerations reveal the similarity of results in this context with these simple toy 
model binaries. In fact, the mean field analysis confirms all the conjectures made earlier based on the binaries. 
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Figure 12. The isothermal cylinder density curve. There is no singularity in contrast to isothermal sphere density. 
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